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Abstract 

We introduce a new paradigm for quantum computing called Ancilla-Driven 
Quantum Computation (ADQC) combines aspects of the quantum circuit [T] 
and the one-way model P| to overcome challenging issues in building large-scale 
quantum computers. Instead of directly manipulating each qubit to perform 
universal quantum logic gates or measurements, ADQC uses a fixed two-qubit 
interaction to couple the memory register of a quantum computer to an ancilla 
qubit. By measuring the ancilla, the measurement-induced back-action on the 
system performs the desired logical operations. 

By demanding that the ancilla-system qubit interaction should lead to uni- 
tary and stepwise deterministic evolution, and that it should be possible to 
standardise the computation, that is, applying all global operations at the be- 
ginning, we are able to place conditions on the interactions that can be used for 
ADQC. We prove there are only two such classes of interactions characterised 
in terms of the non-local part of the interaction operator. This leads to the 
definition of a new entanglement resource called twisted graph states generated 
from non-commuting operators. The ADQC model is formalised in an algebraic 
framework similar to the Measurement Calculus [5]. Furthermore, we present 
the notion of causal flow for twisted graph states, based on the stabiliser for- 
malism, to characterise the determinism. Finally we demonstrate compositional 
embedding between ADQC and both the one-way and circuit models which will 
allow us to transfer recently developed theory and toolkits of measurement-based 
quantum computing directly into ADQC. 

Key words: Models of Quantum Computation, Multi-partite Entanglement. 



1. Introduction 

There are two main paradigms which have driven both the theory and imple- 
mentation of quantum computation; gate-based quantum computing (GBQC) jl], 
and measurement-based quantum computing (MBQC) [2]. Though these two 
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models are computationally equivalent, in practice each has their own advan- 
tages and disadvantages which have major implications for the choice of physical 
system, design, and operation. We introduce a new paradigm called ancilla- 
driven quantum computing which combines features of both models, in order 
to parallelise the architecture of quantum computers, to decrease decoherence 
effects, and simplify their physical implementation and operation. 

GBQC requires, in general, arbitrary networks of dynamic operations which 
in turn complicates the design and characterisation of the entire computer. 
There are many scenarios where it would be desirable to physically separate 
preparation, measurement, and coherent interaction regions to reduce control 
complexity, circuitry congestion, and decoherence due to cross-talk. In contrast, 
MBQC performs computation purely through single-qubit measurement on a 
pre-existing static multi-partite entangled state, distilling and processing non- 
local correlation. However, the generation of the initial highly entangled state, 
incorporation of quantum error correction and fault-tolerance, and individual 
measurement of each qubit are issues in many candidate systems. 

Our new model of ancilla-driven quantum computing (ADQC) is partly in- 
spired by the previous works of Andersson and Oi in ji] , and also Perdrix and 
Jorrand in [5]. Andersson and Oi introduced an efficient method to implement 
any generalised quantum measurement by coupling the system with an ancilla 
qubit. The method, however, assumes arbitrary dynamic coupling operations 
between ancilla and system. Perdrix and Jorrand, on the other hand, describe 
a probabilistic version of MBQC in terms of a Turing machine where one can 
view the read- write head as an ancillary qubit, though this is not the way the 
model is presented. Moreover their approach still requires direct manipulation 
of the memory register and dynamic global measurement operators. 

ADQC attempts to overcome such issues by performing computation where 
the memory register (input data) can only be remotely manipulated through in- 
teraction with a supply of prepared ancillas. In other words, instead of directly 
manipulating data qubits to perform universal quantum logic gates or measure- 
ments, ADQC uses a fixed two-qubit unitary interaction to couple the memory 
register of a quantum computer to an ancilla qubit. By measuring the ancilla, 
the measurement-induced back-action on the system performs the desired log- 
ical operation. Practically, a single fixed unitary interaction coupling the data 
and ancilla qubits greatly simplifies design, construction, and operation of the 
computer since only one particular discrete operation needs to be generated and 
characterised. Furthermore, separating interaction and measurement leads to 
a parallel structure with possibly reduced decoherence [30j . A requisite inter- 
action for ADQC already exists in a variety of physical systems ranging from 
ion micro-traps, neutral atoms, nuclear spin donors in semiconductors, SQUIDs 
and cavity QED which greatly increases the scope for implementation of the 
core ideas. ADQC also naturally benefits from optimisation of the qubit species 
employed for memory and ancilla. Memory qubits can be chosen for long co- 
herence time at the expense of being static and difficult to manipulate directly, 
whilst ancilla qubits may be chosen for high mobility and rapid initialisation 
and measurement, e.g. donor nuclear spins in isotopically pure silicon as mem- 
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ory and electron spins conveyed via charge transport by adiabatic passage as 
ancilla in solid state quantum computing. 

So far we have mentioned only the practical advantages of our proposed 
architecture. The formalisation of the computational model underlying ADQC, 
which is the focus of the current paper, leads to the introduction of a new multi- 
partite entanglement resource. Only recently has it been demonstrated that a 
very restricted class of multi-partite entangled states are useful for universal 
MBQC [SI [7]. However a full characterisation of such states [8] remains an open 
problem which this paper aims to make progress upon. 

The entangled graph states [TD] have emerged as an elegant and power- 
ful quantum resource, especially for measurement-based quantum computation 
(MBQC) [2]. Many important results on their entanglement properties [TT| . 
information flow |13l 114] , implementation |15j , and novel applications in cryp- 
tography [Ml [17], are due to their simple description. The generating operator 
for graph states, called controUed-phase, is a symmetric and commuting oper- 
ator which leads to a simple graphical notation and hence the name for these 
states. Additionally, the elegant result by van den Nest et al. [TI] [TB] shows 
that any stabiliser state is equivalent to a graph state up to local Clifford oper- 
ators. This greatly expands the scope of the results for graph states, and leads 
to a natural extension of the above constructions into stabiliser states, as well 
as allowing a convenient graphical notation for a very general class of states. 
If we consider open graph states, graph states where some nodes (called input 
nodes) are given in arbitrary states (rather than being prepared in a particular 
fixed state which is the case for graph states) much of the theory still follows. 
However open stabiliser states with arbitrary input nodes no longer fulfil the 
pre-requisites of the theorem by van den Nest et. al, and in general they do 
not admit a trivial graphical notation. 

Moreover, not all two-qubit interactions between system and ancillary qubits 
can be used for ADQC. To characterise interactions that will eventually enable 
universal ADQC we demand that the ancilla-system interaction leads to a uni- 
tary and stepwise deterministic evolution of the system qubit, and that it should 
be possible to standardise the computation, that is, applying all global oper- 
ations at the beginning. By doing so we are able to place conditions on the 
possible interactions resulting in two classes of interactions that are necessary 
and sufficient for ADQC. This naturally leads to the definition of a particular 
class of open stabiliser states, called twisted graph states which, despite being 
generated by non-commuting operations, still admits a simple graph representa- 
tion. They form the key ingredient for ADQC. We then show how this new class 
of states can be viewed as open graph states up to somelocal swap operations. 
We also develop an algebraic framework similar to the measurement calculus, 
which is the mathematical framework underlying MBQC computation. This 
makes it possible to derive the standardisation theory for the ADQC patterns 
of computation. As we will see, any ADQC computation requires a classical con- 
trol structure to compensate for the probabilistic nature of the measurement. 
In order to characterise the determinism, we introduce the notion of causal flow 
for twisted graph states based on the stabiliser formalism. Compared to the 
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open graph state, the stabiHser state has a more comphcated and global struc- 
ture. We construct, however, direct translations between ADQC and MBQC 
for a subclass of deterministic patterns with flow which preserve depth, in order 
to prove that ADQC is as parallel as MBQC. We also present the embedding 
between GBQC and ADQC and show how, as for MBQC, a separation in depth 
can be obtained. 



2. Preliminaries 

We briefly review the required concepts from quantum computing. A more 
detailed introduction can be found in [T9|. Let H denote a 2-dimensional com- 
plex vector space, equipped with the standard inner product. Wc pick an or- 
thonormal basis for this space, label the two basis vectors |0) and and 
identify them with the vectors (1, 0)"^ and (0, 1)"^, respectively. A qubit is a unit 
length vector in this space, and so can be expressed as a linear combination of 
the basis states: 



ao|0)-t-ai|l) = 




Here ag^ai are complex amplitudes, and |q;oP + — 1- 

An m-quhit state is a unit vector in the w-fold tensor space TH. (i) ■ ■ ■ ®'H. 
The 2™ basis states of this space are the m-fold tensor products of the states 
|0) and We abbreviate |1) (g) |0) to |1)|0) or |10). With these basis states, 
an 771-qubit state \<j)) is a 2™-dimcnsional complex unit vector 

\4>) = ^ ai\i). 

There exist quantum states that cannot be written as the tensor product of 
other quantum states. Such states are called entangled states, e.g. |00) 4- 

We use ((/)| — 10)* to denote the conjugate transpose of the vector and 
{(p \ ip) — {4)\ ■ 1-0) for the inner product between states |0) and These two 
states are orthogonal if (0 | V') = 0. The norm of |0) is || ||= ^ 4>). 

A quantum state can evolve by a unitary operation or by a measurement. A 
unitary transformation is a linear mapping that preserves the norm of the states. 
If we apply a unitary [/ to a state |0), it evolves to ?7|0). The Pauli operators 
are a well-known set of unitary transformations for quantum computing: 




and the Pauli group on n qubits is generated by Pauli operators. Several other 
unitary transformations that we will use in this paper are the identity 1, the 
phase gate P{c(), of which P{ti /A) and P{tt/2) are a special cases, the Hadamard 
H, the controUed-Z (AZ) and the SWAP operation: 
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The Clifford group on n qubits is generated by the matrices Z, H, P{tt/2) 
and AZ, and is the norniahser of the Pauh group. This set of matrices is not 
universal for quantum computation, but by adding any single-qubit gate not 
in the ChfTord group (such as P(7r/4)), we do get a set that is approximately 
universal for quantum computing |19j . 

The most general measurement allowed by quantum mechanics is specified 
by a family of positive semi-definite operators Ei = M* Mi, 1 < i < k, subject 
to the condition that 'YliiEi — 1. A projective measurement is defined in the 
special case where the operators are projections. Let |0) be an m-qubit 
state and B — {\hi) , . . . , I62'")} an orthonormal basis of the m-qubit space. A 
projective measurement of the state \(j)) in the B basis means that we apply the 
projection operators Pi — \bi){hi\ to The resulting quantum state is \bi) 
with probability Pi = |(0 | An important class of projective measurements 

are Fault measurements, i.e. projections onto eigenstates of Pauli operators. 

So far we have dealt with pure quantum states. A more general represen- 
tation with density matrices also allows us to describe open physical systems, 
where one can prepare a classical stochastic mixture of pure quantum states, 
called mixed quantum states. For a system in a pure state jV'), the density ma- 
trix is just the projection operator Suppose that we only know that a 
system is one of several possible states | V'l ),-.■, | V"?; ) with probabihties pi, . . . ,pk 
respectively. We define the density matrix for such a state to be 



The most general physical operator that acts over density matrices is a com- 
pletely positive trace preserving map (CPTP) £ : B{S)i) — > B{^2) with Kraus 
decomposition 



where the Km ■ S)i S^2, B{:^) is the Banach space of bounded linear operators 
and we require that 



Any unitary operation U can be approximated with a circuit C, using gates 
from a fixed universal set of gates. The size of a circuit is the number of gates and 
its depth is the largest number of gates on any input-output path. Equivalently, 
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the depth is the number of layers that are required for the parallel execution of 
the circuit, where a qubit can be involved in at most one interaction per layer. 
In this paper, we adopt the model according to which at any given time-step, 
a single qubit can be involved in at most one interaction. This differs from 
the concurrency viewpoint, according to which all interactions for commuting 
operations can be done simultaneously. 

2.1. Measurement-based model 

We give a brief introduction to measurement-based quantum computing 
(MBQC) EH [23]. A more detailed description is available in [3l [24 l [26] . 
Our notation follows that of [3]. In MBQC, computations are represented as 
patterns, which are sequences of commands acting on the qubits in the pattern. 
These commands arc of four types: 

1. Ni is a one-qubit preparation command which prepares the auxiliary 
qubit i in state |-|-) = ;;^(|0) + The preparation commands can 
be implicit from the pattern: when not specified, all non-input qubits are 
prepared in the |+) state. 

2. Eij is a two-qubit entanglement command which applies the controlled- 
operation, AZ, to qubits i and j. Note that the AZ operation is symmetric 
so that Eij — Eji. Also, Eij commutes with Ejk and thus the ordering of 
the entanglement commands in not important. 

3. is a one-qubit measurement on qubit i which depends on parame- 
ter a e [0, 27r) called the angle of measurement. M" is the orthogonal 
projection onto the states 

l+a>-^(|0)+e-|l)) 

followed by a trace-out operator, since measurements are destructive. We 
denote the classical outcome of a measurement performed at qubit i by 
TTii G Z2 . We take the specific convention that = if the measurement 
outcome is |-I-q), and that = 1 if the measurement outcome is \~a)- 
Outcomes can be summed together resulting in expressions of the form 

m = rrii 

which are called signals, and where the summation is understood as being 
done modulo 2. The domain of a signal is the set of qubits on which it 
depends (in this example, the domain of m is /). 

4. Xi and Zi are one-qubit Pauli corrections which correspond to the appli- 
cation of the Pauli X and Z matrices, respectively, on qubit i. 
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In order to obtain universality, we have to add a classical control mecha- 
nism, called feed-forward, which allows measurement angles and corrections to 
be dependent on the results of previous measurements [2 [3]. Let m and n be 
signals. Dependent corrections are written as and and dependent mea- 
surements are written as „[M"]™. The meaning of dependencies for corrections 
is straightforward: — — 1 (no correction is applied), while X} = X^ 
and Zl = . In the case of dependent measurements, the measurement angle 
depends on m, n and a as follows: 

„[Mf]"-Aff (1) 

so that, depending on the parity of m and n, one may have to modify the 
angle of measurement a to one of —a, a + tt and —a + tt. These modifications 
correspond to conjugations of measurements under X and Z: 

X™Mf X™ = Mf"^^"" (2) 
Mf = Mf (3) 

and we will therefore refer to them as the X- and Z-actions, or alternatively 
as the X- and ^-dependencies. Since measurements are destructive, the above 
equations simplify to 

MtXr = Mf-i)'"" (4) 
Mf = +"^. (5) 

Note that these two actions commute, since —a + tt = —a — tt up to 27r, and 
hence the order in which one applies them does not matter. 

A pattern is defined by the choice of a finite set V of qubits, two not nec- 
essarily disjoint sets I Q V and O Q V determining the pattern inputs and 
outputs, and a finite sequence of commands acting on V . We require that no 
command depend on an outcome not yet measured, that no command act on 
a qubit already measured, that a qubit be measured if and only if it is not an 
output qubit and that a qubit be prepared if and only if it is not an input qubit. 
This set of rules is known as the definiteness condition. 

A pattern is said to be in standard form if all the preparation commands Ni 
and entanglement operators Eij appear first in its command sequence, followed 
by measurements and finally corrections. A pattern that is not in standard form 
is called a wild pattern. Any wild pattern can be put in its unique standard 
form [5]; this form can reveal implicit parallelism in the computation. The 
procedure of rewriting a pattern in its standard form is called standardisation. 
This can be done by applying the following rewrite rules: 



E,,Xr ^ XrZJ'E,, (6) 

EijZ"'' ^ Z^Eij (7) 

[M^rxf ^ n[M^r+p (8) 

,[Mf]™Z,f => „+p[Mf]'". (9) 
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The rewrite rules also contain the following free commutation rules which tell 
us that if we are dealing with disjoint sets of target qubits, then measurement, 
corrections and entanglement commands commute pairwise [3], so that 

EijA]^ ^ A]f^Eij where A is not an entanglement command (10) 
AkXl^ X™Ak where A is not a correction command (11) 
Ak^™ Z^Ay^ where A is not a correction command (12) 

where k represent the qubits acted upon by command A, and arc distinct from i 
and j. 

Standardisation allows us to graphically present the global operation of a 
pattern. We define an open graph state {G,I,0) to consist of an undirected 
graph G together with two subsets of nodes / and O, called inputs and outputs. 
We write V for the set of vertices in G, E for the set of edges, I'^, and 
for the complements of / and O in 1/ and Eq ■= Hfi j}G-E Eij for the global 
entanglement operator associated with G. Trivially, any standard pattern has a 
unique underlying open graph state, obtained by neglecting measurements and 
correction commands. 

We now give a condition on geometries under which it is possible to synthe- 
size a set of dependent corrections such that the obtained pattern is uniformly 
and strongly deterministic, i.e. all the branches of the computation are equal, 
independently of the angles of the measurements. In what follows, x ^ y de- 
notes that X is adjacent to y in G and iV/c denotes the sequence of preparation 
commands 11 jg/^ 

Definition 1. A flow (/, ^) for a geometry {G,I,0) consists of a map f : 
Qc jc ^ partial order ^ over V such that for all x € O'^ : 

• (^) f{x); 

• (a) xdi f{x); 

• (Hi) for all y ^ f{x), x ^ y . 

The coarsest order ^ for which (/, ^) is a flow is called the dependency order 
induced by / and its depth is called flow depth. 

Theorem 1. Suppose the geometry (G,I,0) has flow f. Then the pattern 

no.. n- (^;'(.) n ^^^^rO^G 

where the product follows the dependency order ^ of f , is uniformly and strongly 
deterministic, and realizes the unitary embedding 

Ugj,o,c. := 2lo^l/2 ( n {+c.^. )Eg. 

If the underlying geometry of a pattern has a flow and its pattern command 
sequence is constructed as given by the flow theorem, we call this pattern a 
pattern with flow. 
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3. Ancilla-Driven Model 

As mentioned in the introduction, in ancilla-driven quantum computing we 
are interested in two essential properties: 

• The only global operation is a fixed two-qubit interaction between any 
ancilla and system qubit. 

• Only ancilla qubits will be measured. 

We introduce the ADQC model within an algebraic framework similar to 
that of the measurement calculus recalled in the previous section. We have a 
set of fixed basic commands described below where the indices i, j, . . . represent 
the qubits on which each of these operations apply. A pattern is a sequence of 
commands defined over a set of qubits in the list V, called computation space, 
where the particular sub-list S represents the system qubits (we may refer to 
them as data or memory register) and the rest A= V\S are the ancilla qubits. 
We define an arbitrary pure single qubit state by 

|+e,^)=cos(f)|0)+e^*sin(f)|l), 
and denote its orthogonal state (the opposite point in the Bloch sphere) with 

he,0)-sin(|)|O)-e''^cos(|)|l), 
where < 6* < tt and < < 27r. 

• Preparation, n]^'^ (a e A) prepares an ancilla qubit in the state It/i). 

• Interaction. Eas (s e S*, a G A) entangle a system qubit and an ancilla 
qubit with interaction operator hZ followed by Hadamard on each qubit: 

KZ:=Hs® Ha hZas 

Note that a SWAP f\Z interaction would be another possible choice. In 
section |4] we will prove that f\Z and SWAP f\Z are the only two possi- 
ble interactions, up to local unitary equivalence, that allow for universal 
ADQC. Choosing one over the other will depend on the natural dynamics 
of the physical system used for an implementation. 

• Ancilla Measurement. Af^'" (a £ A) measures qubit a on plane A G 
{{X, Y), {X, Z), {Y, Z)}, defined by orthogonal projections into: 

- |±(x.y),a) :-|±f,a> if A-(X,y) 

- l±M,a> :=l±a,o) if A=(X,Z) 

- \±iY.z),a) :=|±a,t)ifA=(r,Z) 
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with the convention that \+e.ti>){+e,<t>\a corresponds to the outcome rUa = 
0, while \—e,,p){—e,,j>\a corresponds to = 1. The propagation of depen- 
dent corrections (next command) defines dependent measurement: 

where m,n, . . . are module 2 summation of several measurements out- 
comes, also called signals. The domain of a signal is the set of qubits on 
which it depends]^ 

• Corrections. Xi and Zi {i gV), 1-qubit Pauli operators. As in MBQC, 
to control the non-determinism of the measurement outcomes certain local 
corrections will depend upon previous measurement outcomes. These will 
be written as Cf*, with C° = 1, and C,. 

We write i^y for the associated quantum state space (gj^gvC^. To run a 
pattern, one prepares the system qubits in some given input state ^' G ios, 
while the ancilla qubits are all prepared according to the N commands in fixed 
lip) states. The commands are then executed in sequence, and finally the result 
of the pattern computation is read back from the system qubit^ Similar to the 
one-way model we will consider only patterns satisfying definiteness conditions. 

The main differences between ADQC and MBQC are: (1) The interaction 
operator being AZ instead of AZ, which still belongs to the so-called Clifford 
group, the normaliser of the Pauli group; (2) Only ancilla qubits can be mea- 
sured, that is to say in the terminology of MBQC any ADQC pattern has the 
same number of inputs and outputs which are overlapping (system qubits). 
Apart from universality, which we will prove later most of the theory of mea- 
surement calculus [3] which was developed for the one-way quantum computer 
can be easily adapted to ADQC. For completeness we briefly review this here. 

The first way to combine patterns is by composing them. Two patterns *Pi 
and ^2 may be composed if 5*1 = 5*2. Provided that has as many system 
qubits as by renaming these qubits, one can always make them composable. 
However it is important to emphasise that since the Eij operators are non- 
commuting, their order of appearance in each pattern must be preserved under 
the renaming and composition. The other way of combining patterns is to tensor 
them. Two patterns *Pi and '^2 niay be tensored if Vi O V2 = 0. Again one 
can always meet this condition by renaming qubits in a way that these sets are 
made disjoint. 

3.1. The semantics of patterns 

We present a formal operational semantics for ADQC patterns as a proba- 
bilistic labelled transition system, similar to [3j. Besides quantum states, one 



^Depending on the context we sometimes use the notation m for syntax, i.e., a set of qubits 
(representing a formal sum) and sometimes for semantics, i.e., or 1. 

^Preparation and readout of the system qubits can be performed by using suitable ancilla 
states and measurements. 



10 



needs a classical state recording the outcomes of the successive measurements 
made in a pattern. If we let U stand for the finite set of qubits that are still 
active (i.e. not yet measured) and W stands for the set of qubits that have 
been measured (i.e. they are now just classical bits recording the measurement 
outcomes), it is natural to define the computation state space as 

C := i:u,w^u X Zf. 

In other words, the computation states form a U, VF-indexed family of pairs q, 
r, where g is a quantum state from Sju and F is a map from some W to the 
outcome space Z2. We call this classical component F an outcome map, and 
denote by the empty outcome map in We need further notation. For any 
signal m and classical state F e , such that the domain of m is included in 
W, we take mr to be the value of m given by the outcome map F. That is to 
say, if m = mi, then mr := F(i) where the sum is taken in Z2. Also if 
F S Z^, and 2; G Z2, we define 

T[x/i]{i) = X, T[x/i]U) = T{j) for j ^ i 

which is a map in Z^'^^*\ 

We may now view each of our commands as acting on the state space C: 







9 0|V')i,r 








KZi^q^T 




q,T 




X^g.F 




q,r 

UU{i},W,q,r 
UU{i},W,q,r 




Zf-q,V 

C/,W^U{i},(+A,ar 
f/,VFU{i},(-A,ap 


li9.r[o/i: 
li9.r[iA: 



where ar = (— l)'"^Q; + nr7r. We introduce an additional command called signal 
shifting: 

g,F ^ q,T[T{i) + mr/i] 

It consists in shifting the measurement outcome at i by the amount mr. Note 
that the Z-action leaves measurements globally invariant, in the sense that 
I+c+tt), |— c+tt) = |~a)) l+a)- Thus changing ato a + ir amounts to exchanging 
the outcomes of the measurements, and one has 

Signal shifting allows us to dispose of the Z action of a measurement, sometimes 
resulting in convenient optimisations of standard forms. In the rest of the paper, 
for simplicity, we omit the superscript F on the measurement outcomes. 

The convention is that when one does a measurement the resulting state is 
renormalised and the probabilities are associated with the transition. We do 
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not adhere to this convention here, but instead leave the states unnormaHzed. 
The reason for this choice is that in this way, the probability of reaching a given 
state can be read off its norm, and the overall treatment is simpler. 



3.2. Denotational Semantics 

We now present the denotational semantics of ADQC patterns. If n is the 
number of measurements then the run may follow 2" different branches. Each 
branch is associated with a unique binary string n of length n, representing the 
classical outcomes of the measurements along that branch, and a unique branch 
map (Kraus operator) Kj^ representing the linear transformation from fis to ^35 
along that branch. This map is obtained from the (un-normaliscd) operational 
semantics via the sequence (g^jF^) with 1 < i < m (where m is the total number 
of commands), such that 

gi,ri = q® |+...+),0 

and for all i < m : q,;_i,ri_i — ^ qi,Yi- 

and all measurement commands in the sequence {Ki\ have been replaced by 
appropriate projections corresponding to the outcome index n. 

Definition 2. A pattern *P realizes a map on density matrices p given by p t—^ 
Ks{p)K ■ We write for the map realised by 

It is then easy to prove [5] that each pattern realizes a completely positive trace 
preserving (CPTP) map and if a pattern is strongly deterministic (see section 
2), then it realizes a unitary embedding 3 . Hence the denotational semantics 
of a pattern is a CPTP-map. It is also compositional, as the following theorem 
shows. 

Theorem 2. For two patterns ?Pi and ^2 we have l^i^2\ = I^JsH^Pi] and 

Proof. Recall that two patterns '^2 may be combined by composition 
provided *Pi has as many system qubits as *P2- Suppose this is the case, and 
suppose further that *Pi and respectively realise some CPTP-maps Ti and 
T2. We need to show that the composite pattern ?P2^5i realizes T2T1. Indeed, 
the two diagrams representing branches in *Pi and ^2'- 

•Qsi ^ -Ssi ^ fls2 

y 



can be pasted together, since Si — S2, and 9)Si — ^82- But then it is enough 
to notice 1) that preparation steps p2 in ^2 commute with all actions in *Pi 
since they are appled on disjoint sets of qubits, and 2) that no action taken in 
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^2 depends on the measurements outcomes in It follows that the pasted 
diagram describes the same branches as does the one associated to the composite 
^2^1- A similar argument applies to the case of a tensor combination, and one 
has that ^2 *Pi realizes T2 «) Ti. □ 

3.3. Generating patterns 

In order to prove the universality we present two simple generic patterns 
where only measurements in the {X,Y) plane (M"), Pauli Z measurements 
(M^) and preparations of the ancilla in the state |+) (TV) are sufficient. Note 
that a Pauli Z measurement can be considered as a special case of a measurement 
in the {X, Z) or (Y, Z) plane, with a = 0. 

The following one-parameter family J(a) generates all single-qubit unitary 
operators p7j: 

as any unitary operator U on C'^ can be written: 

U = e'°'J{0)J{l3)J{-f)J{S) 

for some a, (3, 7 and S in M. Recall that the MBQC implementation of the J 
generator is the following pattern: 

3(-a) := X™iMi^-^)^"£;,, (14) 

where s is the system qubit input, a is the ancilla and Ega is AZ operator |27j . 
On the other hand, the following pattern also implements the J gate, but now 
the ancilla qubit a will instead be measured: 

:= i/,Zr"M(^-^)^"£;,„ (15) 

where Hs represents the application of a Hadamard gate on the system qubit. 
We can now manipulate the new pattern to derive a generating pattern for the 
operator J in our model. We want to use Ega as the interaction command: 

a(-a) := i/,Zr"M(^'^)'"£;,„ 
= X^-Mi'^^^^^^HsEsa 

= xr-Afi^'^)^"^,,. (16) 

In addition to this, we only need a generator for a two-qubit unitary such as AZ 
to obtain the full universality. The MBQC pattern for controlled- {Eij with 
i, j € S) is, however, not desirable, as it is an operator between two qubits of the 
system rather than an interaction between system and ancilla qubits. Therefore 
the natural choice instead is to consider the interaction Eas'Eas. It is easy to 
check that this, combined with Pauli Z measurement of the ancilla, will give us 
a simple generating pattern for the two qubit operator AZ: 

A3 :- X'^-M^ Eas'Eas. (17) 
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Any unitary can then be simulated by sequential and parallel compositions 
of the above generating patterns, where the composition simply glues given pat- 
terns over the common system qubits, while preserving the initial orders of the 
commands. We will return to the important issue of how to represent the com- 
posed pattern graphically, but in order to do so, we first have to address the 
important feature of standardization in the ADQC model. That is, the stan- 
dardisation procedure which permits us to rewrite any well defined patterns, e.g. 
obtained from composition, in a standard form where all the preparation com- 
mands are applied first, followed by the entangling, measurement, and finally 
correction commands. 

For simplicity, in the remainder of this paper we will restrict ourselves to 
a special class of patterns, namely, those using only ancillas of degree 1 with 
arbitrary (X, Y) plane measurement and of degree 2 with Pauli Z measurement. 
However, the whole theory developed in this paper can easily be extended to 
the more general setting. 

3.4- Standardisation 

Similar to the one-way model [3] we present a simple calculus of local equa- 
tions by which any general pattern can be put into a standard form where entan- 
glement is done first, then measurements, then corrections. The consequences of 
the existence of such a procedure (called standardisation) are far-reaching and 
is explained in details in [3 . We just recall that since entangling comes first, 
one never has to do "on the fly" entanglements and the rewriting of a pattern 
to standard form reveals parallelism in the pattern computation. In a general 
pattern, one is forced to compute sequentially and to strictly obey the command 
sequence, whereas, after standardization, the dependency structure is relaxed, 
resulting in lower computational depth complexity. It is known that any MBQC 
model can admit a standardisation procedure if and only if the entangling com- 
mand belongs to the normaliser group of the group generated by the correction 
commands [14 . This is the case for our ADQC model. The required rewrite 
rules are 



The rules for propagation of the correction through measurement are the same 
as for MBQC, but with additional rules for the AI^ measurement: 



E^Zl = X-Eij. 



(18) 
(19) 



\M-Y^zi 



(20) 
(21) 




(22) 
(23) 
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We also have the same free commutation rewrite rules: 



EijA^ A]f^Eij where A is not an entanglement command (24) 
A^X^ ^ X™A\s_ where A is not a correction command (25) 
A\s_Z^ => Zj"Ai^ where A is not a correction command (26) 

where k represent the qubits acted upon by command A, and are distinct from 
i and j. 

Recall that the effect of a Z correction on a qubit a simply flips the out- 
come of a measurement to be made on that qubit. Hence we can replace the 
dependencies induced by the Z correction by appropriate operations over the 
measurement outcomes as described below. In what follows, m[n/mi] denotes 
the substitution of rrii with n in m, where m, n are modulo 2 summations of 
several measurement outcomes: 

r,[M^r = K[M:r (27) 

-^mpn _ pn -^m[n+mi/mi] {"iS) 

n[MfrFf = i^f„[p+™./™.][M;]™[^'+'"'/->l (30) 

pmpn _ pnpm[n+mj/mj] j-^^-j 

One can then use the exact same method as for MBQC in order to prove that 
this rewrite system has the desired properties of confluence and termination. 

We emphasise again_that a key difference between ADQC and MBQC is the 
interaction command Eij versus Eij. The explicit inclusion of the additional 
local Hadamard operations, Eij = HiHjEij, is necessary for universality. This 
is due to the fact that no system qubit can be directly measured but instead 
any operation has to be implemented via the ancilla. It is apparent that while 
rotations in the 2-basis can be performed with only E^j = AZij, no basis change 
at the system can be implemented via the ancilla. 

As shown in section 2, for any standard pattern in MBQC we can write 
its underlying open graph state with qubits representing the nodes and Eij 
the edges of the graph. Then, remarkably, only from the geometry of this 
graph we can obtain the dependency structures to guarantee a deterministic 
computation in MBQC. In other words, the simple graph representation for the 
global operation defining pattern allows one to determine dynamic properties 
directly from the static structure. Can we still obtain similar properties for our 
new model? Despite the non-commutativity of Eij the answer is yes. In section 
5 we define the twisted graph state which is the underlying geometry of a given 
ancilla-driven pattern obtained from standardisation and we present how one 
can directly construct the dependency structure from their geometry. 
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4. Characterisation of Interaction 

A central question in the theory of measurement-based quantum computing 
is the characterisation of the universal resources [H [9l [HI [12]. Unlike these 
studies that analyse the computational power of a uniform class of multi-partite 
entangled states, we take a "bottom up" approach and identify basic building 
blocks that can be composed to perform universal ADQC. This requires the 
characterisation of all two-qubit interactions which couple any two ancilla and 
system qubits, while satisfying certain desirable conditions such as stepwise 
determinism, unitarity and standardisation that will lead to universality, in the 
'universal state preparation' sense [TT| [T^]. We obtain a full characterisation 
of universal resources for ADQC for these conditions, while in contrast, such a 
general result is not available for MBQC. 

In this paper we will focus on stepwise determinism, that is a pattern which 
is deterministic after performing each single measurement together with all the 
Pauli corrections depending on the result of that measurement. Other compu- 
tation strategies could be considered that are not stepwise deterministic, such 
as the (finitely) repeated application of the same operation in the scheme |Hlin]. 
Generalising ADQC to adapt to these less restricted strategies remains an open 
problem. 

We introduce the canonical decomposition of two-qubit unitaries |28| . 

Eas = W^. ) Das K®Vs), (32) 

where W'^ , Ws and , Vs are single qubit unitaries and the diagonal global 
operation is 

The vector a = {ax,ay,az} characterizes all non-local properties of Eas- It is 
sufficient to restrict the a's to < a^^y^z < 7r/4, and the set of distinct a, up 
to symmetries, make up the so-called Weyl chamber. The complete character- 
isation of two-qubit interactions that allow for universal ADQC is summarised 
in the following theorem: 

Theorem 3. The global interactions, Eas, between any two system and ancilla 
qubits which enable stepwise deterministic and unitary evolution that also admits 
standardisation procedure for universal ADQC are locally equivalent to (i.e. Das 
is of the form) 

AZ and SWAPAZ. 

Note that, both AZ and SWAP AZ are among the so called "maximally 
entangling" operators, hence the importance of the above theorem is in the only 
if part which proves for the first time the necessity of this particular type of 
building blocks. 

The proof of Theorem 3 is lengthy and we have broken the characterisation 
of the non-local operation Das into several lemmas on unitarity, determinism, 
standardisation and universality. 
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Da, 



Figure 1: The qubit s belonging to the system in state p and the ancilla qubit, a, are 
coupled via the two-qubit unitary Eas- For suitably chosen local unitaries and W^, the 

initial ancilla state V^' |+)a = |+7,(5)a and the measurement basis {W'^ = 
{|+e,c/))a, | — 9,</))a} are such, that the final state of the system is related to the initial state 
by a single qubit unitary on s, = pU^\ conditional on the measurement outcome of 
the ancilla. 



4-.1. Conditions for unitarity 

In order to guarantee a stepwise deterministic computation we require that 
the apphcation of a two-qubit interaction Eas , followed by measurement of the 
ancilla shall result in an effective unitary evolution of a single system qubit, this 
is called unitary condition. In what follows, we absorb the local unitaries W'^ 



and Va in Equation ( 32 1 into the preparation and measurement of the ancilla, 
by choosing appropriate initial state and the measurement basis for ancilla (see 
Figure [T]) : 

K'(l+>a) = \+,,&)a 
W'(l+>a), W^:'(h)a)} - {| + M)a,|-..0>a} 

Lemma 3. The interaction Eas satisfies the unitary condition only if at least 



one tti in Equation (33) is zero. We chose without loss of generality = 0. 
For a^^OLy it is necessary and sufficient that both, the initial state and the 
measurement basis of the ancilla lie in the X — Y plane, i.e. 'y — 9 = tt /2. When 
one a is non-zero, w.l.o.g. ax ^ and ay — 0, then the ancilla parameters must 
obey the following relation: 

sin 9 cos 7 sin (p = cos 9 sin 7 sin S . 

Proof. Depending on the measurement outcome, denoted by + or — , the system 
qubit is transformed according to 



. Ji^, (34) 



where the Kraus operators Kf^ are given by 

Kf = ,(±| Eas \+)a =a {±9 J Ws Das Vs \+^,s)a- (35) 

These Kraus operators shall be proportional to unitaries. 

The local unitaries on the system qubit, Ws and Vs, do not affect whether or 
not the interaction Eas generates a unitary transformation on the system qubit. 
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We will thus focus on the non-local portion Das of Eas for now and reintroduce 
the local unitaries when discussing the standardisation procedure in Lemmas |4] 
and [6] The Kraus operator corresponding to Das alone are 



Kf —a (±61, Das |+7,<5)a 



(36) 



and Das can be diagonalized in the Bell basis [55], 

I'&i) = ^(100) + 111)), |<i>2)- 5^(100) -111)) 
l'i>3) 



^(|01) + |10)), |ci>4) = _(|01)-|10)). 



as 



Dn 



(37) 



with 



7?1 



-a,, 



y^y I u*,2 

ay - a. 



-Olx 



ay 



The completely positive, trace preserving map on p in Equation (34) shall 
have only two possible outcomes, i.e. 



1, = K+Kt^ 



(38) 



The right hand side of the above equation for a general Das and initial ancilla 
state 1+7,5)0 is equal to 

' l + t r 
r* 1-t 



(39) 



where 



t 

r 



sin 2ax sin 2ay cos 7 

sin 7 sin 2az (sin 2ay cos S — i sin 2ax sin S) . 



(40) 
(41) 



Since t — r — must hold, at least one of the a^'s must vanish. 

Let us without loss of generality choose a^ = and hence r = 0. If only 
one a is non-zero then i = is already true. However, if both a.j.,ay ^ 0, 
then the initial ancilla state must lie in the X — Y plane, i.e. 7 = f . These 
two types of interaction are much studied in physics due to their significance 
for coupled spin systems. The case with only one non-zero a is referred to 
as an Ising interaction, and the case with two non-zero a's as a Heisenberg 
interaction |29j . Both cases must have individual Kraus operators that are 
themselves proportional to unitaries. 



KfKf=p±ls 



(42) 
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with probabilities p-^- +p- = 1. This puts requirements on the initial state and 
measurement basis of the ancilla as we discussed next. 

For the Heisenberg interaction with az = 0,ax 0, ay =/= 0, the initial state 
of the ancilla is in the X — Y plane (7 = 7r/2) and Equation (42 1 makes it 
necessary that the measurement basis of the ancilla must lie in the same plane, 
i.e. 6 — Tr/2. The Kraus operators can then be written as, up to insignificant 
phase factors, 



a± —b 



± 



(43) 



with 



S-cj) 



a_ = sm sm ay cos ^ 

, . . S + cf, 

o_ = sm ax cos ay sm — - — 



S-(f> 



I cos ax cos ay sm ■ ^ 



I cos ax sm ay cos ■ ^ 



a-|_ — sm ax sm ay sm 



5-(t) 
2 

5 + 4) 



S-(t> 



(44) 



I cos ax cos ay cos ■ „ 
2 



6+ = sin cos ay cos ^ 



<5 + </) 



I cos Qfa; sm sm ■ ^ 



( 1 ± cos 2ax sin 5 sin (/) ± cos 2ay cos 5 cos < 



The probabilities for obtaining the measurement results + and — are indepen- 
dent of the system p and are given by 

1 

^ 2 

For the Ising interaction with a^ = ay — 0, ax =/= since only one of the oij's 
is non-zero, then any initial state of the ancilla, l+j^sja, will lead to unitary 
Kraus operators for the system qubit provided the measurement basis for the 
ancilla, {\+e.4,)ai \~e.ct>)a}, is chosen such that 

sin 9 cos 7 sin (f> = cos 6 sin 7 sin 5. (45) 

This can be satisfied, for instance, by 9 — j and (j) — 5 for arbitrary initial state 
1+7, (5)a- The Kraus operators for the Ising interaction are given by 

Kf = A±l + i(-l)"±B±X (46) 

where n± are arbitrary integer numbers and the real coefhcients A± , B± are 

cos Q!t 



A± = 
B± = 



V2 



sm a 



± cos 7 cos 6* ± sin 7 sin 9 cos{S — (f)) , 
-\/l+ COS7COS0 ± sin 7 sin cos ((5 + 0), 



where the angle 9 and the phase <j) of the measurement basis has to satisfy 
Equation M5|. The probabilities for the two measurement outcomes are 



1 



(1 ± sin 9 sin 7 cos S cos </> ± cos 2ax (cos 6* cos 7 -|- sin sin 7 sin (5 sin </>)) . 



This concludes the proof of lemma. □ 
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4-2. Conditions for correctable branching 

As the measurement results are random, the unitary operation that is apphed 
to the state is either [/+ or U~ where — / y/p± are the normahsed 
Kraus operators. To generate a deterministic evolution, for example with 
we require that the other unitary, U~ , can be corrected with an additional Pauli 
correction Pg (up to an unimportant relative phase A), 

U-^e^^P,U+. (47) 

General Pauli corrections are of the form Ps{a,b,c) — aXs + bYg + cZg with 
a, 6, c e R and + 6^ + = 1 and have the properties 

P^p\ p2^1, Tr[P]=0. 

Pauli corrections are the standard choice in measurement-based computation. 
Their structure is well understood and essential in ensuring an overall deter- 
ministic computation. Reuniting the branches with a Pauli correction after the 
measurement places no direct constraint on the interaction. Das- Instead this 
can be achieved by choosing the initial state and the measurement basis of the 
ancilla appropriately, as specified below. 



In the Heisenberg case the correctable Kraus operators ( 43 1 must fulfil the 



relation a*^a_ + = 0, i.e. the two non-zero a's must relate to each other 
as 

cos 2ax tan (5 = tan (/)Cos 2Q;y. (48) 

It can be seen that all values of ax and ay can be covered independently by 
choosing the phases S and </> for the ancilla appropriately. 

However, if one of the interaction parameters is chosen as ay — -k /A the above 
relation implies 5 — Q and the probability for each branch becomes p± = |. 
It will be discussed in Lemma |4] that this requirement is indeed necessary to 
allow for the standardisation of a computational pattern. The Kraus operators 



Equation ( 43 1 then have the coefficients 



a_ = -j= ysm ax cos + * cos ax sm 2 ' ^ 

o_ = I sm ax sm — + t cos ax cos — 
v2 \ 2 2 



(49) 



It can be seen that these parameters imply that the Pauli correction between 
the two Kraus operators is Y (up to a phase), = iYs K~ . 

In the Ising case, the Kraus operators are a sum of identity and Pauli-X only. 



see Equation (46) and the only non-trivial correction available is X. To fulfil 
— Xg—p^, one has to choose the ancilla's initial state and measurement 
basis such that 

yl+A_ + (-l)"++"-S+S_ = 0. (50) 

It can be seen that one n has to be even and the other one odd in oder to satisfy 
this condition. The probabilities for the two measurement outcomes need not 
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be balanced and are given as p± = A^. + B\ . The interaction strength ax can 
then be manipulated by choosing the ancilla parameters, 



, / 1 — (cos7Cos6' + sin7sin0cos(i5 — 0))^ 

tan^a^r^W ^ ^ (51) 

y 1 — (cos 7 COS 9 — sin 7 sin 9 cos((5 + (j))) 

4-3. Conditions for standardisation 

A computation can be standardised, see 3.4 when a general Pauli correction 
transforms, under the interaction Das, into a tensor product between ancilla and 
system qubit, where the new correction on the system qubit is again a general 
Pauli operation. By symmetry of Das the same must be true for the ancilla, i.e. 
we require 

Das la(^Ps{a,b,c) D[^^Ta®Qs, (52) 

where Ta and Qg are either general Pauli operations or the identity operation. 
If this relation is valid the correction Pg can be commuted through future in- 
teractions and shifted to the very end of the computation. 

Lemma 4. For the Heisenberg interaction there are two possible solutions for 
pairs of unitaries Das{ax, oty, 0) and Pauli corrections P{a, b, c), up to relabeling 
of the a's, that obey the standardisation relation Equation (52) 

1. Fixed Heisenberg interaction: ay = 7r/4 and a^ = 7r/4. 

Here Pauli corrections P(a, &, 0) and P(0,0, 1) transform according to 
la (E> {aXs + bYs) ^ {aYa - bXa) ^ Zs, a^ + b^^ 1, 

la(E>Zs ^ Za^ls- 

2. General Heisenberg interaction: ay — 7r/4 and a^ ^ 7r/4. 
Pauli corrections P(1,0,0) transform according to 

l„(g)X, ^ Ya®Zs. 

Proof. For Das{oixTay,Q) the general Pauli correction P{a,b,c) transforms 
according to 



PsDas ~ la® (a COS 2ayXs + &COS laxYg + ccos 2ax cos 2ayZs) 

(53) 



-sin2Q!j. Xa ® {~bZs + ccos 2ayYs) 



+ sin 2ay Ya (g) {aZs — ccos 2axXs) 
+ sin 2ax sin 2ay Za® c\s- 



To determine which interactions fulfil Equation ( 52 ) we take the partial traces 
of either the ancilla or the system. Reading from the right hand side of Equa- 
tion (52 1 the traces must be, up to global phases, 

Tr[Ta] • = I Q iff Ta is a general Pauh operation, ^^^^ 

T ■ Tr[Q 1 = i ^"^^ iff = 1 C55') 
a v^s\ -y g iff is a general Pauli operation. ^ ' 
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Figure 2: Composition of gates on a single qubit s. In this example, each line represents 
the same effective computation. Attempting to implement K~^{ipi) results, with probability 
p_, in the actual operation K~{ipi) = P K'^ (ipi) where a correction P appears (first line). 
The correction can be shifted past local unitaries W and V (second line) and past subsequent 
ancilla-driven operations K^(ip2) (third line), see the text for details. 



The partial trace over the ancilla in the Heisenberg case, Equation (53 1, is 



Ti-a[Das'ia^PsDas] — 2 ■ {a COS 2ayX g + b COS 2axYs + c COS COS 2ay Z g) . (56) 

For a^-, cXy this expression can never be equal to 2Qs where Qs is a general 
Pauli or the identity for the system qubit. However, the right hand side can 
vanish (and Ta must hence be a general Pauli) when a — and ax = 7r/4 or 
6 = and ay = 7r/4 or ax = ay = 7r/4. Thus, for the Heisenberg interaction 
one of the parameters has to be fixed to 7r/4 to enable full standardisation of the 
computation and we choose, without loss of generality, ay = tt/A. In the case 
that 6 = 0, i.e. for corrections P{a,0,c), the right hand side of Equation ( [53| 
becomes 

Ya (8i {aZs — ccosaxXg) + csin2aa; Za(^l, (57) 



which implies either c = to match condition Equation (54). When both 



a = 7r/4 then the right hand side of Equation (53 1 becomes 



i-bXa + aVa)® Zs + cZa^l, (58) 

and the correction must thus be either P{a, b, 0) or P(0, 0, 1). □ 

An example for Case 1 is the SWAP AZ interaction whereas Case 2 allows 
a range of parameters < ax < 7r/4; however, only one type of correction, X, 
obeys the standardisation condition. This becomes an issue when composing 
several ancilla-driven operations after one another. 

Lemma 5. Case 2 in Lemma^is not compositional. 

Proof. In order to implement arbitrarily complex computations it is necessary 
to repeatedly apply ancilla-driven unitary transformations to a single system 
qubit, see Figure [2] ^ 

The full ancilla-system interaction Eas plus measurement of the ancilla gen- 
erates the effective operator sequence {±\a Eas \+)a — Wg Kf{ip) Vg on the sys- 



tem qubit, where Ws and Vg see Equation (32 1 are the local unitaries on the 
system qubit and ip denotes all ancilla parameters. Due to the probabilistic na- 
ture of the measurement, when attempting to implement the operation K'^{Lpi), 
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the actual outcome could be K~{ipi) oc P K~^{(fi) with a correction P, pictured 
in the first line, Figure |2] The branching behaviour discussed in the previous 
subsection requires that P — Y for the Heisenberg interaction and P — X 
for the Ising interaction. The local unitaries affect this original correction and 
transform it into a possibly different general Pauli P', as shown in the second 
line, Figure [2j 

Constructing an arbitrary computation requires composition of subsequent 
ancilla-driven operations, K^(ip2) etc. When commuting the correction P' with 
the K~^{ip2) the result must be a new Kraus operation K^{(f'2), implemented 
by a (possibly) different set of ancilla parameters (p'2, and a new general Pauli 
correction P" after the new Kraus operation, as depicted in the third line of 
Figure [2] In the two cases of interest here, cases 2 and 4, Pauli X is the only 
allowed correction, i.e. P' — X. Moreover, in order for the new correction P" 
to be shifted further, it has to be identical to the branching relation P, or 1. 
The pictorial shifting conditions imply the following mathematical relations 

VWP^ e'^ P' VW, 

^"(^^):P-e-' P"4^, (59) 



P = e P or e 1, 



where A, A', A" are arbitrary global phases. 

In Case 2 of Lemma |4] this set of equations cannot be fulfilled, as there is 
no p air o f ancilla parameters ip2 and for which the Kraus operators Equa- 



tion (43l with coefficients Equation (49l obey K+{ip2)Y = e*^ X K'^{ip2) or 
K~^{(p2) Y = e^^ K^{ip'2). This implies that the possible ancilla-driven compu- 
tations using a general Heisenberg interaction are extremely limited - only one 
ancilla interaction per system qubit can be tolerated. To enable composition 
with the Heisenberg interaction forces the interaction parameter to take the 
special value = 7r/4. □ 

Similarly, we characterise the Ising interaction cases. 

Lemma 6. For the Ising interaction there are two possible solutions for pairs 
of unitaries Dasiax,0,0) and Pauli corrections P{a,b,c), up to relabeling of the 



a's, that obey the standardisation relation Equation (52) 



3. Fixed Ising interaction: — 7r/4. 

Here Pauli corrections P(1,0,0) and P(0,fe, c) transform according to 

la^ibYs + cZ,) ^ Xa'Si {bZs - cYs), a2 + 52^1. 

4-. General Ising interaction: 7^ 7r/4. 

Pauli corrections P(1,0,0) transform according to 

Proof. The argument is similar to the proof of the previous lemma. □ 
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Figure 3: Any sequence, U{ipi, ip2, ■ ■ ■) = W U{ipi, ip2, ■ ■ ■)V , of K operations and VW 
operations commutes or anti-commutes with X. 



Case 4 again allows only one type of correction, an X correction. We prove 
next that this requirement restricts any obtained ancilla-driven unitaries on a 
single system qubit to be on a fixed plane and thus fails universality. 

4- 4- Universality 

The general Ising interaction, Case 4, fulfils the set of equations [59] as the 
Pauli X correction commutes with the actual Kraus operators 

K+{ip2)X = XK+{^2) 



for ~ j4+ 1 + i (—1)"+ _B_|_ X, given in Equation (46 1, where the product 
of the local unitaries has to obey VWX = e^^XVW, i.e. VW is either of the 
form (up to global phases) VW = al + ibX or of the form VW = aY + bZ with 
some coefficients a, 6 e IR and + ~ 1. However, this leads to the following 
lemma. 

Lemma 7. Case 4 ^''^ Lemma^is not universal. 

Proof. Any sequence of consecutive application of K{ip) and VW (see Figure|3| 

U{^i,^2,---) = ■■.WK{^3)V WK{^2)V WK{ipi)V 

contains a kernel, U {ipi, ip2i ■ . ■) — W U{ifi, ip2, ■ ■ ■)y ■, that is again of the form 

U{ipi,ip2, • • ■ ) = a(^i, ■ • . ) 1 + ibi^f 1,^2, .■■)X 
or U{ifi,ip2,---) = a{ipi,ip2,-.-)Y + b{ipi,ip2,...)Z. 

The coefficient functions a(i^i, (^2, • ■ • ) and 6(</'i, f2, ■ ■ ■) with a? + b'^ — \ will 
depend on the choice of all ancilla parameters Lpi^if2, ■ ■ ■ in the sequence. It is 
easy to show that both classes of operations, for arbitrary parameters, can move 
a given input state only in one plane of the Bloch sphere, parallel to the Y — Z 
plane. Finally, we note that the additional first and last local unitary operation, 
V and W , which individually can take arbitrary form as long as their product 
VW (anti-)commutes with X, are fixed. That means that once W an V are 
given, the plane where JJ operates is tilted into an arbitrary, yet fixed direction 
in which the state can evolve. Since it is impossible to move out of that plane 
it is implied that the U{(pi, (p2, ■ ■ ■) for general Ising interaction. Case 4, does 
not allow universal single qubit preparation. □ 

The final step for the proof of Theorem [3] is given below. 
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Lemma 8. The remaining cases 1 and 3 are universal. 

Proof. It suffices to identify suitable local operations that ensure universality. 
The hZ interaction is an example for Case 3 and choosing the local unitaries 
Vs and Ws appropriately, for instance: 

Eas = Ha® Hg t\Zas 

- HaP{Tr/2)aHa®HsP{Tr/2)sHs (7r/4, 0, 0) Ha®H, 

makes this model universal as we proved in Subsection |3.3[ The SWAP AZ 
interaction is an example for Case 1 and it is easy to verify that ADQC model 
with SWAP f\Z is indeed the same as the one-way model, as the position of 
ancilla and system qubits are swapped, and hence universal. □ 

Note that the local unitaries on the ancilla are irrelevant as they just redefine 
initial state and measurement basis of the ancilla. Moreover, as said before, for 
the Case 3 the local unitaries on the system have to enable a "change of basis" , 
this would rule out the class of diagonal unitaries. However, the full charac- 
terisation of local unitaries that makes each cases 1 and 3 universal remains 
an interesting open problem. For the rest of this paper we will focus on the 
Eas = Ha ® Hg l\Zas interaction which defines a new class of entangled states. 

5. Twisted Graph States 

The main issue with the Eij operators is the fact that they are non-commuting, 
therefore after standardisation their order will be important. Another impor- 
tant property of an ancilla-driven pattern is that system qubits interact only 
with ancilla qubits. Therefore we introduce a multipartite entangled state as 
a graph over ancilla and system qubits and edges with extra condition to 
address the above mentioned requirements. In Section |6.2| we show how this 
new class of states can be viewed as open graph states up to some local swap 
operators. This is the reason behind the chosen name. 

Definition 9. An open twisted graph state (G, S, A,C) consists of a bipartite 
graph G over disjoint sets of qubits S and A, called systems and ancillas, such 
that the maximum degree of any ancilla node is 2, together with an edge labelling 
C. The labels define a partial ordering over edges where the order is total and 
strict on any edges that share a common vertex, i.e. it defines an edge colouring. 

The corresponding quantum state, denoted as \Eq) , is obtained by preparing 
qubits in S in given arbitrary states and all the qubits in A in the \+) state, and 
then applying Eij over corresponding qubits according to the partial order of C 
(see Figure^. 

One may think of an open twisted graph state as the beginning of the def- 
inition of an ancilla-driven pattern, where one has already decided how many 
qubits will be used {V = S U A) and how they will be entangled: 

Eg :== Il{t,j}eEEij- 
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Figure 4: An open twisted graph state where the system qubits are the white circles that will 
not be measured and the rest are ancilla qubits. The edges are E interactions and edge labels 
denote the partial order. 



To complete the definition of the pattern it remains to decide which angles will 
be used to measure ancilla qubits and, most importantly if one is interested 
in determinism, which dependent corrections will be applied. Conversely, any 
ancilla-driven pattern has a unique underlying open twisted graph state obtained 
by forgetting measurements and corrections where the colour of the edges is 
given by the most partial order of the Eij commands which respects the non- 
commuting order. Recall that two Eij commands commute if and only if they 
act on disjoint set of qubits. Note that the depth of this partial order is the 
true depth of the preparation of the state. 

Different ordered colourings over the same graph structure might lead to dif- 
ferent twisted graph states and consequently different patterns of computation. 
We leave as an open question whether one can find a more relaxed definition 
that can still uniquely define the entangled state corresponding to an ancilla- 
driven pattern. On the other hand our restrained definition will allow us to 
derive the dependency structure of the measurements directly from the order of 
the colouring, this is the topic of the next subsection. 

5.1. Dependency structure 

We will use the graph stabiliser formalism [lOl [19] to construct a determin- 
istic pattern. Recall that for any open graph state \Eg) defined over a graph 
G with vertices V, we have the following set of equations for all the non-input 
qubits: 

X^U,eGi^)Zd\EG)) = \EG) 

where G{i) is the set of neighbour vertices of i in G. The above Pauli operators 
are called the stabiliser operators of \tfj). 

Similarly we define the stabiliser operators of a given twisted graph state 
\Eg) defined over a graph G with vertices V = S U A. We will use the following 
notation as well. Define S{a) for a G A to be the attached system qubit s G S 
with the smallest edge label and S'{a) to be the other one if it exists; N{s) for 
s £ 5 to be the set of ancilla qubits connected to the system qubit s; and finally 
Gs(a) to be the sub-graph with edges between S{a) and Ng(^a.)- 

Consider first a simple case where Eg = £-05(0)^!^^ • Then the stabiliser 
has the form 

ZaXsia){\EG)) = \Eg) (60) 
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The above equation is due to 



and for another simple case of \Eq) = EaS'(a)EaS{a)Nl^'^ we have 

XaXsia){\EG)) = \Eg) (61) 

again due to 

XaXs{a){EaS'{a)EaS(a)Nl'^'^) = EaS' {a)EaS{a)XaZ s{a) Z S{a)Nl'^'' 

~ EaS'(a)EaS{a)Nl'^\ 

In order to generahse the above cases the following rewrite rules for Pauli 
commutations will be used: 



Equation 18 EijXf = XjZ^Eij, transforms the X operation on the sys- 
tem qubit to the next immediate ancilla qubit, introducing a Z operation 
at the system qubit. 



Equation 19 EijZf = XfEij, replaces the Z operation on the system 



qubit with an X operation. 

Unlike the stabiliser of the graph state which has a local structure, in the case of 
a twisted graph state the stabiliser of a may affect the whole of the graph. This 
action is, however, recursive and can be defined using the collection of several 
local actions^ Define the label of an ancilla node to be the same as the label 
of the edge EaS{a)- Consider an ancilla qubit a and those qubits in N{S{a)) 
with label greater than label of a (see Figure |5]). We can assume, without loss 
of generality, that all edges connected to S{a) have labels 1 to n with 1 being 
the label of a. This is due to the fact that the stabiliser of a has an effect 
only on those qubits in N{S{a)) where their edge interaction are after the edge 
interaction of a and S{a) hence having a greater edge label. 



3 n 

Figure 5: The generic case for the local study of stabiliser at ancilla a. 



Definition 10. Given a twisted graph state (G, S, A,C), a local stabiliser on 
the ancilla qubit a is defined as following: 
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• Consider vertices in Gg(^a) with labels greater than label of a and relabel 
them from I to n according to C ordering, with 1 being the label of a. 

• Add X on vertices in N(S{a)) with even label and degree 1. 

• Add Z on vertices in N{S{a)) with even label and degree 2. 

• Add X on the system qubit if n is odd otherwise add Z. 

We denote the above set of Pauli operators with Pi{S{a)) which acts on a subset 
of qubits in N{S{a)). 

Define I{a) to be the set of degree-two ancilla qubits in N{S{a)) with even 
label that get a Z Pauli operator in the definition of the local stabiliser of a. 
The stabiliser of a will have the same local effect as defined above over S{a') for 
all the a' G I{a) and the same structure repeats for vertices in I{a'). Therefore 
we define recursively such qubits 

r{a) = {a'\3n : a' e /"(«)} 

where I^{a) — I{a) and /"''"^(a) = Ua'e/"(a) -^('*')- We can now present a 
recursive definition for the stabiliser of a twisted graph state as the product of 
a collection of local stabilisers. 

Definition 11. Given a twisted graph state (G,S,A,C), the stabiliser on the 
ancilla qubit a is defined as follows: 

P{a) = ZaY{^,(,i.i^^)Pi{S{a')) if a has degree one 
P{a) = XaYl^,^j,(^^-^Pi{S{a')) otherwise 

It is a straightforward but cumbersome computation to show the correctness of 
the above definition and we omit the details of the proof. 

We can now adapt the notion of flow and generalised flow of graph states 
[T^ [T3] for twisted graph states to derive a sufficient condition for determinism. 
The key idea is exactly the same as in the MBQC case based on the follow- 
ing simple observation. We could make a measurement M^^'^'^'''" "determinis- 
tic" (corrected) if it could be pre-composed by an anachronical Z^' correction 
(i.e. conditioned on the outcome of a measurement which hasn't happened yet). 
This unphysical scenario is a useful starting point for our proof. 

The flow construction guarantees that a deterministic pattern with anachronical 
corrections _ 

V - ifaeA {+iX,Y)Ja Eg 

= UaeA M(^'^)'"zr" Eg 

can be transformed into a runnable pattern, where all dependencies will respect 
the proper causal ordering. The key observation which allows us to transform 
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this into a runnable pattern is that the flow construction defines a stabiliser 
Pf{a) which when composed with the anachronical correction, forms an operator 
which commutes with the measurement, and thus the pattern can be brought 
into runnable order. 

For simplicity in the rest of the paper we consider only patterns where degree- 
two vertices are measured with Pauli Z and degree-one vertices are measured in 
the {X,Y) plane, we use the generic term M^"'"" for both cases. This class of 
patterns are large enough to introduce a universal ADQC model as they include 
the generating pattern introduced in Section [3. 3 1 However the definition of flow 
and determinism can be extended to the more general case as well. 

Definition 12. An open twisted graph state (G, S, A, C) has causal flow ij there 
exists a partial order > over V such that for all a & A and all vertices a' G P{o) 
we have a < a' except for those a' that will he measured with Pauli Z . 

Theorem 4. Suppose the open twisted graph state (G, S', A,C) has a causal flow 
with the partial order >. Define: 

C{a) ~ P{a)Za for all degree-one ancilla a 
G(a) — P(a)Xa for all degree-two ancilla a 

then the pattern: 

Vg,c. ■■= UaeA Ciar- M>>"» Eg 

where the product follows the dependency order >, is runnable, uniformly and 
strongly deterministic. 

Proof. The proof is based on the following equations first consider the {X, Y) 
measurement case for degree-one ancillas: 

i+alaiEc) = M-Z^^^{Eg) 

= M-Z^^-P{ar^{EG) 
= Ciar-^M^iEG) 

Similarly for the Z measurement over degree-two ancillas we have: 

{OUEg) = M^XJ^^^IEg) 

= A/fXr»P(ar"(^G) 

= C{ar^M^{EG) 

Hence we can write: 

UaeA{^a,aa\aEG^UaeA C^"" ^'^a"'"" 

The left hand side is clearly a uniformly and strongly deterministic pattern. 
The right hand side pattern is runnable as the introduced corrections follow the 
partial order > except for the Z correction introduced over degree-two ancillas. 
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However one can ignore them since these qubits will be measured with Pauli Z 
and we have Mf Z™ = Mf . This finishes the proof. □ 

It is interesting to note that the flow deflnition for a graph state was based 
on the geometry of the underlying graph, whereas in a twisted graph state it 
is based on the edge colouring order. Indeed, as mentioned before, different 
edge colourings lead to different twisted graph states and hence different flow 
constructions. Roughly speaking, the edge colouring plays the role of geometry 
for the twisted graph states. 



6. Compositional Embedding 

One of the main foci in constructing direct translations between models is 
to study parallelism as the more parallel the computation, the more robust it 
is against decoherence. Recently the advantage of MBQC over GBQC in terms 
of depth complexity has been demonstrated where a logarithmic separation 
was shown [30]. We will prove a similar result for ADQC and we present a 
transformation between ADQC and MBQC that preserves depth. 



6.1. GBQC and ADQC 

The question of translating GBQC circuits into MBQC patterns and vice 
versa has been addressed before in [30j and it can be directly adapted for ADQC 
as well. In fact the universality proof of ADQC already presents a method of 
translation of a circuit into ADQC: (I) Rewrite the given circuit in terms of the 
universal gates set of J{a) and AZ; (II) Replace each gate with its correspond- 
ing ADQC pattern (equations 16 and 17l; (III) Perform the standardisation 
procedure. 

The above construction cannot be used in reverse as the edge colouring 
order might lead to a circuit with an acausal loop. However it is possible to 
keep all the auxiliary qubits to avoid creating loops in the resulting circuit. The 
scheme is simply based on the well-known method of coherently implementing 
a measurement [30]. It is also easy to prove, in a similar way as in [30], that 
the translation from a GBQC circuit into an ADQC pattern will never increase 
the depth as the edge colouring number of the obtained twisted graph state 
will be upper-bounded by the depth of the original circuit. More importantly, 
we present an example where the depth decreases exponentially. Consider the 
ladder structure of the circuit in Figure [6] which has depth n. This circuit, 
through the introduced construction, will be translated into a pattern with the 
twisted graph state shown in Figure |6] which has constant depth 4. This is due 
to the fact that the preparation depth for any ADQC pattern by deflnition, is 
upper bounded by the edge colouring of the graph, which in this case is equal to 
4. On the other hand the Pauli X measurements on the right most qubits, break 
the dependency chain between measurement i.e. an measurement depends 
only on the result of some Pauli measurements. Hence one can flnd a flow for 
this twisted graph state where all the Pauli measurements are in the flrst layer 
of the corresponding partial order and all the non-Pauli measurements are in 
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the final layer and computation depth which is equal to the flow depth is 2. 
It is easy to verify that the same logarithmic depth separation result between 
GBQC and MBQC obtained for the parity function [30] is also valid for the case 
of GBQC and ADQC. 




Figure 6: A ladder structure circuit with the binary l\Z and the unitary P{cti) gates, together 
with the corresponding open twisted graph state obtained via gate by gate translation. The 
edge labels present the edge colouring and the vertex labels are the measurement angles, where 
X and Z stands for Pauli measurement. 



6.2. MBQC and ADQC 

As mentioned before there exists a compositional embedding from MBQC 
patterns with flow into GBQC and vice versa, and together with the construction 
of the last subsection one can obtain an embedding between ADQC and MBQC 
for patterns with flow. However in view of parallelism, it is interesting to flnd 
such an embedding directly by presenting the correspondence between twisted 
open graph states and open graph states. 

The following equation relates the two resources but it is only valid for 
degree-one ancilla qubits 

EasNl+^ = SWAP,, EasNl+^ . (62) 

This equation and the next one are in fact the reason behind the chosen name 
for this class of states as one can recover a graph state from them by applying 
the appropriate sequence of twist (SWAP) operators. In order to handle the 



31 



degree-two ancilla qubits we will use the following pattern equations 



A3 - Xl'^'^M^Eas'Eas 

= Xl^^M^ HaHs'Eas'Eas 
— X^'^M^ Hs'Eas'Eas 

= X^'X^'M^Mi'Ebs'Eas'Eas (63) 
In the new pattern for A3 both instances of Eas can be replaced using Equation 



^ X^^'X^' M^M^ E,s'Ess'Eas (64) 

Therefore we can replace any pattern over a given twisted graph state where 
degree-two vertices are measured with Pauli Z, into a pattern over a graph state 
obtained through the above manipulations of the Eas edges. 

The other direction of translation, i.e. from graph sates to twisted graph 
states can be obtained from the following general equation. 

Ess' = Hs Ess' 

= Hs' Hs X"'"'M^ Eas'Eas 

= X^^'M^Es'b' Xl^'M^Esb XT^M^Eas'Eas 

= M^M^M^ Es'b'Eas'EsbEas (65) 

Not that the resulting twisted graph state might not be unique as it depends on 
the order of re-write rules application. However for the special class of graph 
state with flow T3], which are still universal for MBQC, we present a unique 
translation that is also more efficient, i.e. uses fewer ancilla qubits. Consider a 



graph state (G, /, O) with flow (/, ^), see Appendix 2.1 for definitions. Denote 
the orbit of / on i g / with F.^ = {v £ V\3n : f"{i) = v} and let Oi be the 
unique element in O H Fi. Due to the property of the flow |13j . every v € V 
belongs to a unique Fi so we can deflne the rank of v to be the unique value 
r{v) such that (i) = v . The corresponding twisted graph state Gt{S, A) is 
constructed by the following steps, where the system qubits in Gt are the input 
qubits in G. 

1. For any i £ S, add r(oi) many degree-one ancilla qubits to A and attach 
them to i, and set the edge labels 1, • • • , r{oi). 

2. For any non-flow edges in G between a £ Fi and b G Fj add a degree-two 
ancilla qubit g to ^ with edges Fig and Ejg. Set the edge label of Eiq and 
Ejq to be r(a) -I- 1 and r(b) + 1 (if a and b have the same rank make one 
of the edge label bigger). Finally adjust the labelling of the added edges 
in the flrst step in accordance of these new labels. 
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It is a straightforward but cumbersome computation to show that the above 
translation from MBQC into ADQC is depth preserving so we omit the details 
of the proof. We leave as an open question whether a general depth preserving 
translation between MBQC and ADQC can be constructed. 

7. Discussion 

ADQC presents significant advantages over GBQC for particular physical 
implementations. By isolating the system memory from measurement and state 
preparation, the physical layout of a quantum computer can be optimised. Po- 
tentially decoherent read-out mechanisms can be located away from the mem- 
ory. Since only a fixed two-qubit unitary gate has to be implemented between 
memory and ancilla qubits, this simplifies considerably the construction, char- 
acterisation, control and operation of the computer. Control line clutter can 
be reduced, due to relaxing the need to implement single qubit rotations and 
measurement on the memory, simplifying the architecture of the computer as 
well as minimising the possibility of cross-talk. The choice of physical qubits 
for memory and ancilla can also be optimised. Memory qubits can be chosen 
for long coherence times at the expense of being static, whilst ancilla qubit can 
be chosen for mobility, and ease of initialisation and measurement. 

Natural candidate systems consist of an array of static qubits addressed by 
flying qubits or a "read/write" head. For example, an optical lattice of neutral 
atoms can be addressed by a separately controlled atom [3TJ [52]. Since it is 
diflicult to individually address with lasers a single site of a fully filled optical 
lattice, the read- write head would interact with selected sites and can be used as 
the ancilla. A similar idea can be applied to neutral atoms trapped in dipole trap 
arrays [35] controlled by an atom in an optical tweezer [33]. Similarly in [35] . 
ions in an array of micro traps would be manipulated by a single ion read- write 
head. Alternatively in [36], the role of the memory is played by the quantised 
electromagnetic field in an array of cavities whilst Rydberg atoms traversing 
the cavities act as ancillas (similar ideas are contained in [37]). All of the above 
schemes possess a natural t\Z (dispersive) interaction between the ancilla and 
system qubits. For universality, either an additional Hadamard operation on the 
system should be incorporated into this interaction (e.g. through the pushing 
laser in [35]), or an effective SWAP operation be found in conjunction with 
the ControUed-Z. The latter could be achieved through cold collisions between 
ancilla and system PHI [35] . 

A system which may prove particularly amenable to our model is discussed 
in [30]. Here, the nuclear spin of a single dopant atom in isotopically pure 
silicon plays the role of a memory qubit which can be controUably coupled via 
the hyperfine interaction to an electron spin which acts as an ancilla qubit. 
Nuclear spins can be very well isolated from the environment, as well as the 
state preparation and measurement areas. Electrons can be rapidly transported 
around the computer using charge transport via adiabatic passage (CTAP), 
this avoids the issue of swapping nuclear spin states, as in the original Kane 
proposal [4P, which can lead to a reduction in fault tolerance. An issue here 
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is making sure that the interaction between electron and nuclear spins is of the 
correct form as to allow conditional unitary dynamics. An interaction of the 
form l\Z followed by SWAP gate can be achieved by using the method presented 
in |42j with only the Heisenberg interaction between ancilla and system and local 
operations on the ancilla itself. 

A natural interaction which is also suitable for ADQC is the XY-Hamiltonian 
which be easily turned into the SWAP f\Z gate g3j (equivalent to the ISWAP 
in the aforementioned reference). In |44j . a natural XY-interaction between 
nuclear spins, as in the Kane proposal, is mediated by a 2D electron gas in 
the Quantum Hall regime. The XY-interaction also naturally occurs between 
quantum dots coupled by a cavity |l5] or superconducting qubits [15J |37] . Here, 
the memory and ancilla qubits are of the same species. 
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